Dynamic van der Waals Theory of two-phase fluids in heat flow 
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We present a dynamic van der Waals theory. It is useful to study phase separation when the 
temperature varies in space. We show that if heat flow is applied to liquid suspending a gas droplet 
at zero gravity, a convective flow occurs such that the temperature gradient within the droplet 
nearly vanishes. As the heat flux is increased, the droplet becomes attached to the heated wall that 
is wetted by liquid in equilibrium. In one case corresponding to partial wetting by gas, an apparent 
contact angle can be defined. In the other case with larger heat fiux, the droplet completely wets 
the heated wall expelling liquid. 

PACS numbers: 05.70.Ln, 64.70. Fx, 44.35. +c, 47.20. Dr 



O 



I 

o 



> 

ON 

o 

o 
vn 
o 



o 



In usual theories of phase transitions, the fluctuations 
of the temperature T are assumed to be small and are 
neglected. However, there can be situations in which 
phase transitions occur in inhomogeneous T. For exam- 
ple, wetting properties near the gas-liquid critical point 
are very sensitive to applied heat flux [1,2] and boiling 
processes remain largely unexplored [3,4]. To treat such 
problems we propose to start with coarse-grained entropy 
and energy, from which the temperature is defined as a 
functional of the order parameter and the energy density. 

Let us consider one-component fluids, where the or- 
der parameter is the number density n. In the van der 
Waals theory the energy density e and the entropy s 
per particle are given by e = dnkBT/2 — evori^ and 
s = (d/2)lnr -I- ln(l/won — 1) + const., respectively, in 
terms of T and n [5] . Here vq and e represent the molec- 
ular volume and the magnitude of the attractive poten- 
tial, respectively, and d is the space dimensionality. To 
describe situations in the presence of interfaces, we gener- 
alize this description by including gradient contributions 
to the total entropy and the total internal energy as 
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The gradient terms represent a decrease of the entropy 
and an increase of the energy due to inhomogeneity of n. 
Here the internal energy density is written as 



e + K\V7i\^/2, 



(3) 



including the gradient part. We assume that the coeffi- 
cients C and K are simply constants and that s in Eq.(l) 
depends on n and e as 

s = fee ln[(e/n -f evQn)'^/'^{l/von - 1)] + const. (4) 



This relation follows from the van der Waals theory. Our 
scheme starts with Eqs.(l), (2), and (4), where the temer- 
ature is not yet introduced. We then define the local 
temperature T = T(n, e) by 



(5) 



where n is fixed in the derivatives. This definition of T 
is analogous to that in a micro-canonical ensemble. We 
also define the local chemical potential fi (per particle) 
including the gradient contributions by 
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where e in Eq.(3) is fixed, /i = —T(d{ns)/dn)e, and 
M ^CT + K. 



(7) 



Maximization of S under a fixed total particle number 
J dm and a fixed total energy £ leads to the equilibrium 
conditions T =const. and ft =const. As first derived by 
van der Waals [5,6], the equilibrium interface density pro- 
file n = n{x) is determined by /i(n, T) — Md^n/dx"^ — /icx 
with the surface tension 7 = M dx{dn/dx)^ , where 
/icx is the chemical potential on the coexistence curve [7] . 

Next we set up the dynamic equations [8]. Hereafter 
gravity will be neglected. The mass density p — mn (m 
being the particle mass) obeys dp/dt — —V • (pv), where 
V is the velocity field. The momentum density pv obeys 



d_ 

dt 



pv ~ —V • {pvv) — V • (n — a). 



(8) 



The H = {Hy} is the reversible stress tensor invariant 
with respect to time reversal v —v. The ct = {cr.y } is 
the dissipative stress tensor of the form = r]{\7iVj + 
jVi) -f (C — 2ri/d)5ijV ■ V with — d/dxi in terms 
of the viscosities ?/ and Q. The (total) energy density 
ex = e -|- pv^ /2 including the kinetic part is governed by 
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— CT = -V • [e^v + (il - ct) • ■«] + V • (AVT), (9) 

where A is the thermal conductivity. We construct the 
stress tensor 11^ such that the entropy production rate 
within the fluid assumes the usual form [8], 

= / dri[V • AVT + ^ a,, V,t;,]. (10) 

If there is no heat flow from outside, the above quan- 
tity becomes non-negative-definite. To obtain Eq.(lO) 
we need to require the general relation, 

Using the van der Waals pressure p = nkBT/{l — vqu) — 
evon? we then find 

Hij = {p + Pi)Sij - MViTiVjTi, (12) 

where pi = M{nV'^n + |Vn| V2) + TnVn • \/{M/T). 

As an example, we consider a bubble with radius R 
in heat flow. Thermocapillary bubble migration has long 
been observed especially in space [9-11]. Theoretically, 
however, the bubble motion has been examined in the ab- 
sence of first-order phase transition at the interface such 
as in the case of air bubbles in silicone oil [9]. In such 
cases the temerature dependence of the surface tension 7 
causes a bubble velocity of order vy = R"/i{dT/dz)oc/'ri 
in a given temperature gradient {dT/dz)oo (Marangoni 
effect), where 71 = —dj/dT > and 77 is the viscosity 
in liquid. In one-comonent fluids the effect is drastically 
altered due to latent heat generation or absorption at the 
interface [12]. We point out that a convective velocity of 
order vq = X{dT /dz) 00/ p'T As inside the bubble is suffi- 
cient to carry the applied heat flux, where p' is the mass 
density in gas, As is the entropy difference between gas 
and liquid per unit mass, and A is the thermal conduc- 
tivity in liquid. Note that the ratio 

Ml = vy/vo = Rjip'TAs/rjX (13) 

is a very large number of order R/a with a being a mi- 
croscopic length. With phase change the bubble velocity 
is of order vq and is very slow [12]. Furthermore, the 
temperature inside the bubble should be nearly homoge- 
neous due to the latent heat convection. This is a natural 
consequence in one-component fluids without contami- 
nation [1,12], because the pressure should become nearly 
homogeneous outside the bubble and the fluid near the 
interface should be close to the coexistence curve in the 
p-T phase diagram. 

In our simulation we set K = and then the inter- 
face thickness is i = {Cvo/^kBY^^ far from the critical 
point (since {p, — Pc^)/k-QT is of order 1 for n ^ v^^). 
Space is measured in units of t and the system length is 



L = 202£. Time is measured in units of = ml'^/rivo. 
Here we assume constant viscosities with r] = (. The 
strength of dissipation is represented by the normalized 
viscosity rj* = v^rij l^fnu.. We set rf = 0.1^/^. On all the 
boundary walls we assume [13,14] 

n = riw + ivfV ■ Vn, (14) 

where v is the outward normal unit vector at the sur- 
face. We set riw = 5nc/2 and — 2^/21, for which 
the wall is wetted by liquid in equilibrium in our simu- 
lation. We prepared an initial equilibrium state with a 
circular gas bubble with radius 50 placed at the center 
of the cell. The temperature was at T 0.875Tc, where 
the density is O.STric in gas and 1.74nc in liquid. The 
bottom boundary [y = 0) was then increased by a con- 
stant AT, while the top boundary {y = L) was held at 
the initial temperature Tt = 0.875Tc. (Here we use "bot- 
tom" and "top" though no gravity is assumed.) The side 
boundaries (x = and L) are insulating [v ■ VT = 0). 
Furthermore, the density-dependence of the thermal con- 
ductivity is crucial away from criticality [15], so we set 
A = {vQrj/kBm)n with the thermal diffusivity Dt being 
of order 1 (in units of £^/to = voTj/m). Then A in liquid 
is larger than that in gas by 5 times. 

(i) First, for very small AT, the bubble position is 
shifted toward the bottom without shape changes. Fig.l 
shows such a steady state at AT = 0.00675Tc taken at 
t = 3.4 X 10^ after application of heat flux. Here the 
migration is stopped by the repulsion from the walls due 
to the wetting condition in Eq.(14). The velocity field 
V is almost parallel to the interface, but there is a small 
velocity component through the interface. Thus, at the 
interface, first-order phase transition takes place and la- 
tent heat is generated or absorbed. As demonstrated in 
Fig. 2, latent heat transport is so efficient that the tem- 
perature gradient almost vanishes within the bubble. 

(n) Second, for the case AT = 0.054Tc, Fig.3 shows 
droplet migration toward the bottom. Fig. 4 displays the 
velocity field in the steady state taken at t = 8 x 10^. 
The velocity component through the interface is more 
increased than in Fig.l, again resulting in a flat temper- 
ature inside the droplet as in Fig. 5. As can be seen in 
Figs. 4 and 5, a thin liquid layer is sandwiched between 
the bottom boundary and the droplet. The thickness of 
this layer is so thin that we can define an apparent con- 
tact angle ^eff) which is a decreasing function of AT. In 
accord with these results, Garrabos et al. observed that 
gas spreads on a heated wall initially wetted by liquid 
and exhibits an apparent contact angle [1]. 

(iii) Third, wo increased AT to 0.10125Tc, where the 
bottom temperature is 0.97Tc. As shown in Fig. 6, phase 
separation occurred to produce a gas domain at the bot- 
tom, into which the suspended droplet was finally ab- 
sorbed. In the steady state in Fig. 7, the bottom is cov- 
ered by a gas layer or is completely wetted by gas [16], 



but the interface with the bulk hquid is stiU curved and 
a small velocity field is induced. As can be seen in Fig. 8 
taken at t = 3 x 10^, the temperature gradient within 
the gas layer becomes steeper than that in the liquid. It 
is well-known that a gas film on a heated wall suppresses 
heat transfer to the bulk liquid and, as it thickens, film 
boiling is induced in gravity [3]. 

In the steady states the velocity is of order vq ~ 
DtAT/TL in terms of Dt = X/pCp from Cp ~ As. 
As a result, the Reynolds numbcr(^ Pr~^RAT/LT) is 
very small, where the Prandtl number Pr = voTj/rnDT is 
here taken to be of order 1. In fact, in Figs.l, 4, and 7, 
the variance y/Jv^ is 0.12, 1.9, and 0.63, respectively, in 
units of 10~^£/tQ, where £/to = DtIPtI. In our simula- 
tion the viscous dissipation is very strong even in tran- 
sient states, while in the theory in Ref. [11] the inertia 
effect (vapor recoil) is important in bubble motion on a 
heated wall. This difference stems from the fact that our 
droplet size is very small if t. is microscopic. 

In the diffuse-interface methods for crystal growth [17] 
and thermo-capillary flow [18,19], the interface width I 
may be conveniently treated as a free parameter. This 
is because the phase field equations yield the moving- 
interface problems with appropriate interface boundary 
conditions as long as I <C -R(= domain size). Also in 
fluids we may take £ as a free parameter in the range 
£ <^ R as far as we treat the steady state profiles (where 
the Marangoni effect is suppressed and the Stokes ap- 
proximation is valid). However, to describe rapid fluid 
motions, the normalized viscosity r?*(oc r]/£) should be 
made very small if £ is enlarged. 

In summary, we have developed a scheme to study 
phase transition dynamics with inhomogeneous temper- 
ature under mass and heat fluxes. It can serve as a 
phase field model in fluids, though we should examine 
its sharp- interface limit [19]. For one-component fluids 
we mention wetting dynamics with evaporation and con- 
densation [2,13], as exemplified here, and boiling or con- 
vection in gravity [4,5]. These effects can in principle be 
studied in the scheme of thermocapillary hydrodynamics 
supplemented with interface conditions. However, hydro- 
dynamic calculations accounting for latent heat transport 
have rarely been performed. 
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AJfT^ = 0.00675 
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FIG. 1. Suspended gas droplet in a steady state in heat flow with AT = 0.0067571;. The arrows indicate the velocity. 
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ATfTc = 0,054 




t=1000 t=7000 t= 80000 

FIG. 3. Migration of a gas droplet for AT — 0.054Tc. It apparently wets the bottom partially in the steady state (left), 
where a thin liquid layer is still sandwiched. 





FIG. 5. The temperature in the nearly steady state in Fig. 4. 
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ATfTc = 0.10125 




t= 300 t= 800 t= 30000 

FIG. 6. Time evolution for AT = 0.10125rc. Phase separation occurs at the bottom. The bottom is completely wetted by 
jas in the steady state. 



AT/T^ =0.10125 
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FIG. 7. The velocity field around the gas region at t = 3 X 10" in Fig.6. 



AT/Tc = 0.10125 




FIG. 8. The temperature in the nearly steady state in Fig. 7. 
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